Henry's law, surface tension, and surface adsorption in dilute binary mixtures 
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Equilibrium properties of dilute binary fluid mixtures are studied in two-phase states on the basis 
of a Helmholtz free energy including the gradient free energy. The solute partitioning between 
gas and liquid (Henry's law) and the surface tension change A7 are discussed. A derivation of the 
Gibbs law A7 = —TV is given with F being the surface adsorption. Calculated quantities include the 
derivatives dT c /dX and dp c /dX of the critical temperature and pressure with respect to the solute 
molar fraction X and the temperature-derivative (dj /dT)cx,-p of the surface tension at fixed pressure 
p on the coexistence surface. Here (dy /<9T) CX , P can be both positive and negative, depending on 
the solute molecular size and the solute-solvent interaction, and diverges on the azeptropic line. 
Near the solvent critical point, it is proportional to (dp c /dX) / Kkt, where Kkt is the Krichevskii 
parameter. Explicit expressions are given for all these quantities in the van der Waals model. 



I. INTRODUCTION 

Many problems in physics and engineering involve di- 
lute solutions. In one-phase states, the critical behavior 
of dilute fluid mixtures have been studied extensively-i^i 3 -, 
where crossover occurs from pure-fluid behavior to 
binary-mixture behavior on approaching the critical line. 
In two-phase states, it has been of great interest how a 
solute is partitioned between gas and liquid and how it 
is adsorbed in or repelled from the interface region^&i. 

In the dilute limit, the solute-solute interaction may 
be neglected for nonelectrolytes. Nevertheless, the two- 
phase behavior is still highly nontrivial, depending sen- 
sitively on the detail of the solute-solvent interaction. In 
particular, the surface tension change A7 due to a solute 
is related to the excess solute adsorption 8 . To under- 
stand such effects, we will present a simple Ginzburg- 
Landau theory of dilute mixtures including the gradient 
free energy 3 -. As is well-known^, van der Waals originally 
constructed such a theory for pure fluids to describe a 
gas-liquid interface and to calculate the surface tension 
7. For binary mixtures it is moreover possible to calcu- 
late the solute density profile around an interface, which 
should satisfy the Gibbs adsorption law. In this approach 
solute partitioning between the two phases may be exam- 
ined systematically. 

In fluid hydrodynamics involving a gas-liquid interface, 
it is crucial how the surface tension varies on the sur- 
face as a function of ambient temperature, concentration, 
and pressure, since its variation induces a Marangoni 



,10.11 



However, the present author is not aware of any 



fundamental theory on the surface variation of 7 in fluid 
mixtures in nonequilibrium. Hence we will also calculate 
the surface-tension derivative (dj/dT) cx ^ p with respect 
to the temperature T at fixed pressure p in two-phase 
coexistence^. 

In Section II, we will present a Ginzburg-Landau model 
to calculate how the coexistence surface and the crit- 
ical line are formed with addition of the second com- 
ponent. Mean-field critical behavior of dilute mix- 
tures will be discussed, where the so-called Krichevskii 
parameter— i^i^ i 13 i 14 ' 15 will be of crucial relevance. On 



the basis of the Gibbs adsorption law to be derived in 
Appendix A, general expressions for the surface tension 
variations on the coexistence surface will be given. In 
Section III, use will be made of the van der Waals free 
energy of dilute mixtures^^ supplemented with the gra- 
dient free energy. It will give explicit expressions for all 
the physical quantities discussed in Section II, in terms of 
two dimcnsionless parameters characterizing the solute- 
solvent interaction. In Appendix B, correlation- function 
expressions for thermodynamic derivatives including that 
of the Krichevskii parameter will be giver 3, 1 



II. THEORETICAL BACKGROUND 
A. Ginzburg-Landau theory 

This paper treats dilute nonelectrolyte binary mixtures 
with short-range interactions undergoing the gas-liquid 
transition. The number densities of the two compo- 
nents are written as n\ and rii with 712 <C n\, which 
are coarse-grained variables changing smoothly in space. 
A Ginzburg-Landau free theory is used to describe two- 
phase coexistence. A number of authors calculated the 
surface tension of mixtures by combining an equation of 
state and the gradient theor y 18 ' 19 i 20 . 

Hereafter the Boltzmann constant will be set equal to 
unity. The free energy functional F = F{m, 712} depends 
on n\ and n-i as 



F 



dr 



f- 



E 



Dij\Jrii ■ \/rij 



(2.1) 



The first term / = f(n-\_,n%,T) in the brackets is the 
Helmholtz free energy density dependent on the densities 
and the temperature T. The gradient terms are needed 
to account for a free-energy increase due to density inho- 
mogeneity. The coefficients -Dn, -D12 = -D21, and D22 are 
assumed to be constants independent of the densities. In 
the dilute case n 2 the following form is assumed: 

/ = /o(n x ,T) + Tn 2 [ln(n 2 A 3 ) - 1 + ip{n u T)]. (2.2) 
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Here the van der Waals attractive interactions among 
the molecules of the species 2 (ex rzj) are neglected. 
The fo(ni,T) is the Helmholtz free energy density of 
the one-component (pure) fluid of the species 1 and 
A2 = h(2n/m2T) 1 ^ 2 (h being the Planck constant) is the 
de Broglie length of the species 2. The term Tn 2 f arises 
from the solute-solvent interaction, where ip — ip(ni,T) 
is independent of rii (see the next section for its van der 
Waals expression). 

For the free energy density / in Eq.(2.2) the chemical 
potentials of the two components (without the gradient 
contributions) are expressed as 



/'1 

/'2 



df_ 
dm 
df_ 

dn 2 



(io(ni,T) + Tn 2 p'(ni,T), (2.3) 
Tln(n 2 Xl)+Tip(n u T), (2.4) 



where Uq = dfo/drii is the chemical potential of the pure 
fluid and ip' — dip/dni in fjti. Note that /i2 tends to —00 
logarithmically in the low density limit n 2 — * 0. The 
pressure p = nifii + n 2 p: 2 — f is expressed as 



p = ni/io - fo + Tn 2 (l + nup'), 



(2.5) 



where the last term is the solute correction. For the free 
energy functional F in Eq.(2.1) the generalized chemical 
potentials including the gradient contributions read 



the second line. In the second line n 2 is expressed in 
terms of ri\ and V 2 ni. It follows n 2 — C exp[— f(n\, T)] 
in the homogeneous bulk region. In our theory expan- 
sions up to first order in £ or f 2 are performed. On the 
other hand, Leung and Griffiths 1 used another parameter 
Clg = + A exp(/ii/T — n 2 /T)] in order to describe 
the overall thermodynamics of binary mixtures along the 
critical line (0 < X < 1), where Aq is an appropriate 
constant. 

Equilibrium states may be characterized in terms of 
the field variables, T and £, (instead of T and the average 
solute density). As a functional of ri\ parameterized by 
T and £, the grand potential is defined as 



n = F 



dr(flini + jj, 2 n 2 ). 



(2.10) 



In the dilute case jl 2 may be removed with the aid of 
Eqs.(2.4) and (2.6), leading to 



n= dr 



T 

fo - Ai«i + ^-£>n|Vni| 2 



Tn? 



(2.11) 



where the gradient terms proportional to D12 cancel to 
vanish and n 2 depends on n% as in the second line of 
Eq.(2.9). Here f2 = f2({ni}, T, £) is minimized in equilib- 
rium as a functional of n\ In fact Sfl/Sni = holds from 
Sn 2 — ~n2[f' — Di2\7 2 ]Sni at fixed T and 



% 3=1,2 



D l3 V 2 nj (£ = 1,2), (2.6) 



which are homogeneous in space in equilibrium. The 
usual chemical potentials [i\ and fi2 deviate from /ti and 
fl 2 in the interface region. Originally, van der Waals set 
up the following interface equation for pure fluids 9 , 



fio(n,T) - TD n n" = /i° x (T), 



(2.7) 



where /i° x is the chemical potential on the coexistence 
curve p = p®y.{T) of the pure fluid. The density n(z) 
changes along the z axis and n" = d 2 n/dz 2 . Our equa- 
tions in Eq.(2.6) lead to the van der Waals interface equa- 
tion (2.7) for m = n(z) and n 2 = 0. 

In this paper a small parameter £ is defined as 



C = A ; 



-3JVT 



(2.8) 



which has the dimension of density. The solute density 
n.2 is expressed as as 



B. Two-phase coexistence 

Let a planar interface separate gas and liquid regions. 
The bulk densities of the two components far from the 
interface are written as nn, nx g , 1121, and ri2g- The sub- 
scripts £ and g stand for liquid and gas, respectively. This 
paper treats the dilute regime, 



n 2g < nig, n 2 £ <C n u , 



(2.12) 



in the two phases. Hereafter thermodynamic relations in 
this case are given. For a noncondensable g solute, 
another typical situation is given by n\ g < ri2g nu far 
below the solvent criticality. 

As a reference state, we consider the two-phase state 
of the pure fluid composed of the first component at the 
same temperature T below T c q, where n% is equal to nio 
in liquid and n g o in gas. The chemical potential and 
pressure in the pure fluid are written as /i° x and p° x , 
respectively. With addition of solute, Eq.(2.9) yields the 
bulk solute densities, 



n 2 = Cexp 
= Cexp 



- ip(ni,T) + D 12 V 2 ni + D 22 V 2 n 2 

- tp(m,T) + -Di 2 V 2 ni . (2.9) 



The fugacity of solute f 2 = exp(/i 2 /7 1 ) = \ 2 ( is usu- 
ally used to represent the degree of solute doping. The 
term proportional to X7 2 n 2 in the first line is omitted in 



n 2l = Ce~ 



n 2g = (e~ 



(2.13) 



where (p a = (p(n a o,T) with a standing for t or g. Since 
the pressure is given by the common value p = p cx in the 
two phases, Eq.(2.5) yields the coexisting solvent densi- 
ties nia (a = £ or g) as 



1 — K Ta 

n a o 



Spcx ~ Tn 2a (l + n a0 ip'(n a0 )) 



, (2.14) 
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where Kx a is the isothermal compressibility Kt — 
{dn/dp)T/n of the pure fluid for a — i or g and Sp cx = 
Pcx — Pc X is the deviation of the coexisting pressure. Fur- 
thermore, Eq.(2.3) yields 



Mlc 



/x° x = (8p cx - Tn 2a )/n a0 (2.15) 



for the deviation of the solvent chemical potential in two- 
phase coexistence. This holds both for a = I and <?, so 
Sp cx (l/n m - l/rigo) - T(n 2e /n e o - n 2g /n gQ ) = 0. Thus, 



S Pcx = TAX/ Av, 
Su lcx = TAn 2 /An, 



(2.16) 
(2.17) 



where An = neo — n g o and Av = 1 /n g o — 1 /neo are the 
differences of the density and the volume (per particle) 
between gas and liquid in the pure fluid, respectively, 
(taken to be positive). The differences of the solute den- 
sity and molar fraction are written as 



An 2 = n 2g - n 2e = (< 



AX 



n 2g _ n2£ 
n g0 n eo 



-<Pb - e - 




( 'erf a 


e -<pt 


\ n g0 





C, (2.19) 



which are both proportional to C from Eq.(2.13). 

For infinitesimal variations of fii, p 2 , T, and p, the 
Gibbs-Duhem relation generally holds in the form, 



-XdA - sdT + vdp, 



(2.20) 



where A — [i 2 — fi± is the chemical potential difference, 
s is the entropy per particle, and v = l/(ni + n 2 ) is the 
volume per particle. In particular, for variations on the 
coexistence surface in the p-T-A space, we obtain 



AXdA = -AsdT + Avdp, 



(2.21) 



where As and Av may be taken as the entropy difference 
of the pure fluid. Here dA = T(~ x d( in the dilute case, 
so in the mixture case £ > we have 



dT 



= -T 



.AX 
CA? 



AX 
= T 



(2.22) 



where {d ■ ■ ■ /d ■ ■ -)c X ,p and (d ■ ■ ■ /d ■ ■ -) C x,t arc the 
derivatives on the coexistence surface at fixed p and T, 
respectively, and the right hand sides of Eq.(2.22) are in- 
dependent of ( since X cx Obviously, Sp cx in Eq.(2.16) 
follows from integration of (dp/dC) cx ,T in Eq.(2.22) with 
respect to C from the reference pure fluid state at fixed T. 
To derive <5/ii cx in Eq.(2.17) we integrate Eq.(2.20) with 
respect to ( at fixed T to obtain Eq.(2.15). Here note the 
relation XdA ^ T J Q C d(X/( = TX, where X/C is 
independent of C- In the same manner, the temperature 
change ST CX = T cx (p 7 ^) — T® x (p) at fixed p (below the 
critical pressure p c o) on the coexistence surface reads 



8T CX = -TAX /As, 
which is proportional to 



(2.23) 



It is convenient to introduce the partition coefficient 
of solute K. as the ratio of the solute molar fraction in 
gas X g = n 2g /n g o and that in liquid Xi — n 2 i/nicr£. 
Equation (2.13) gives 



X g n m 
K, = — = exp(( y a f - tp g ). 



(2.24) 



Then AX = (1 - K^Xg = (JC - l)X e . The azeotropic 
line on the coexistence surface is determined by K, = 1, 
on which the two phases have the same composition. If 
the gas region is dilute, X g is nearly equal to the partial 
pressure of the second component divided by the total 
pressure in the gas region. Near the critical point K, — > 
1. When the gas phase is dilute, Henry's constant H is 
usually defined as 



n = P2g /x e 



(2.25) 



withp2g(— Tn 2g ) being the partial pressure of the solute, 
Here TL = p g IC, where p g is the total gas pressure. To 
analyze data near the critical point Sengers et alA used 
another definition of Henry's constant, 



k H = f 2 /Xt 



(2.26) 



where f 2 is the solute fugacity. In our notation we obtain 
k H = KngoXle*" from Eqs.(2.8) and (2.9). 



C. Surface tension and surface adsorption 

The surface tension 7 in binary mixtures may be calcu- 
lated from Eq.(2.1). It has been calculated in the gradi- 
ent theory in fair agreement with experimental data over 
a wide temperature rang o 18 i 19 ' 20 i 21 . However, our result 
cannot be used in the asymptotic critical region. 

In Appendix A, the deviation A7 = 7 — 70 will be 
calculated, where 70 is the surface tension in the pure 
fluid. For small £ it follows the Gibbs relation^, 



A 7 = -IT. 



(2.27) 



Here T is the excess adsorption of the solute on the in- 
terface expressed as 



dz 



An 2 

n 2 (z) - n 2l + -^{n{z) - neo) 



where An 2 



n 2a 



A? 

n 2 i and An 



neo 



>y0 



(2.28) 



and the 



integrand is nonvanishing far from the interface. 

The physical meaning of T is as follows. For a finite sys- 
tem with length L much longer than the interface width, 
the interface position z = z- m may be determined with 
the aid of the Gibbs construction, 



neo + (L- 
Then T is expressed as 



dzn(z). 



(2.29) 



J dz[n 2 {z) - n 2e ] + j dz[n 2 (z) - n 2g ], (2.30) 
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where the first (second) term represents the excess ad- 
sorption in the liquid (gas) region. The integrands here 
tend to far from the interface, so we may push the lower 
bound in the first integral to — oo and the upper bound 
in the second integral to oo for a macroscopic system. 
The Gibbs relation (2.27) has been used frequently for 
surfactants added in water-air and water-oil systems 22 , 
which induce a dramatic decrease of 7 even at extremely 
low bulk densities. If salt is added, A7 contains an elec- 
trostatic contribution also 2 ^. 

The surface tension 7 = j(T, C) of mixtures is defined 
on the coexistence surface p = p cx (T, () . Since A7 oc 
C in Eq.(2.27), use of Eq.(2.22) gives the temperature 
derivative of j(T, £) at fixed p in the form, 



97 
df 



drfy 
dT 



As 
AX' 



r. 



(2.31) 



It is important that the second term is independent of £ 
as well as the first term. In the azeotropic case AX = 0, 
the second term in the right hand side tends to ±00. 



D. Mean-field critical behavior 

1. Landau expansion 

The mean-field critical behavior of dilute binary mix- 
tures will then be examined near the critical point of the 
pure fluid (solvent criticality) . The critical temperature, 
pressure, and density at the solvent criticality are written 
as T c o, PcOj and n c o, respectively, in the pure fluid. The 
order parameter is the solvent density deviation, 



ip = m- n c0 . 



(2.32) 



Here T — T c q and ip are assumed to be small. The Landau 
expansion of /o(«i, T) is of the form, 



A 



fo = hc{T) + ii 0c {T)^+'-f(T-T cQ )^ + ^\ (2.33) 

where f Qc (T) = f (n c0) T) and n 0c (T) = fi (n c0 ,T) are 
the free energy density and the chemical potential at the 
critical density, respectively. The Gibbs-Duhem relation 
for one-component fluids yields 



Voc(T) = n cQ - (s c0 - n cQ p' cy .){T - T c0 ) 



(2.34) 



where /i c o is the critical chemical potential, s c q is the 
critical entropy, and p' cx — (<9p/<9T) cx is the derivative 
of p with respect to T along the coexistence line at the 
solvent criticality. Use has been made of the relation 
(dp/dT) n = p' cx near the solvent criticality^. 

A small amount of the second component is then added 
as a solute. Near the solvent criticality, we expand the 
solute density n 2 in Eq.(2.9) as 



n-2 



G 



C 



[C + d1> + + -^1 (1 - D 12 V 2 # (2.35) 



Here we may set T = T c o, since the term —Tn 2 is already 
a small perturbation in the grand potential (2.11). The 
coefficients Co, C\, C2, and C3 are obtained from the 
expansion of e~ v as 



Co = 
C 3 = 



C 1 = - i p' c C , C 2 = fr? - rf)C , 



Mrt-rt'-ip^Co, 



(2.36) 



where ip' c , tp", and ip" 1 are the derivatives dip/drii, 
d 2 ip/drii, and d^ip/dnf at the solvent criticality. The 
critical solute density and molar fraction read 



n 2c = CCo, X c = (C /n c0 . 



(2.37) 



Equilibrium is obtained by minimization of the grand 
potential f2 in Eq.(2.11), which is the integral of the den- 
sity a/ = fo — fiifii — Tn2 plus the gradient term. Here 
fix should be expressed in terms of the macroscopically 
given pressure p (not treated as a fluctuating variable), 
temperature T, and £. From the expression for /ii in 
Eq.(2.3) some calculations give 



fJ-i = McO 



P-Pco 
n c0 



s c0 (T - T c 



1 



-n 2c , 



(2.38) 



in the bulk regions. This relation also follows from in- 
tegration of the Gibbs-Duhem relation (2.20) for mix- 
tures. Note that \X\ is equal to the right hand side of 
Eq.(2.38) in the whole space. The Landau expansion of 
uj = fo — Ai n i ~ Tn 2 is now of the form, 

w = -po(T, - n~ l h^ + ^(T - T C )V 2 + • ■ • , (2.39) 

where Pq{T,Q — p(n c07 T,() is the pressure in Eq.(2.5) 
at n = n c o and h has the meaning of the ordering field. 
Use of Eqs.(2.34) and (2.38) gives 

h=p- Pc0 - p' cx ( T ~ T co) - T c0 (C - ncoCiX, (2.40) 



In the third term of Eq.(2.39) T c = T c[ 
critical temperature with the shift, 

AT C = (T c0 C 2 /Ao. 



AT r is the 



(2.41) 



Since h = at the criticality T = T c and p — p c , the 
critical pressure shift Ap c — p c — p c o is calculated as 



A Pc = p' cx AT c + T c0 {C - n c0 Ci)C. 



(2.42) 



to first order in £. Since AT C , Ap c , and X c are all linear 
in C, the derivatives of T c and p c along the critical line 
are given by dT c /dX = AT C /X C and dpJdX = Ap c /X c . 
The critical line is characterized by X = X c (() in 
Eq.(2.37), leading to 



dT c „ C 2 

dX ~ nc ° Tc0 Ao~C~o' 



(2.43) 
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In addition, from the third order term (oc C3) in the 
expansion of n 2 in Eq.(2.35), there arises a small shift of 
the critical solvent density as 



n\ c - n c0 = (T c0 C 3 /B . 



(2.45) 



If we expand u) = /o — /*ini — Tn 2 up to the quartic term 
and rewrite it in powers of n\ — ni c , the third order term 
should vanish. However, this critical density shift does 
not affect the shifts of T c and p c to first order in £. Also 
the coefficient of the gradient term in £1 is changed from 
Du to 

D' n =Du -(D 12 d. (2.46) 
This correction is irrelevant in the dilute limit. 

2. Krichevskii parameter and concentration fluctuations 

In the lit er at ur 1 3 1 1 4 1 1 5 ; use has been made of the 
thermodynamic derivative (dp/dX) n T with n = n\ + n 2 
and X = 112/ n to analyze the critical behavior in dilute 
mixtures 13 . From Eq.(2.5) it is equal to Tm(l + n\tf') — 
nf/g in our approximation. It is known to tend to a 
well-defined limit, called the Krichevskii parameter, as 
Q — > at the solvent criticality. In terms of Cq and C\ in 
Eq.(2.41), it is expressed as 

K Kr = (%) = T c0 n c0 (l - Had/Co) (2.47) 



dX 



nT 



From Eqs.(2.43) and (2.44) it follows the well-known 
relation^ 3 -, 



K Kl = 



dp c 
dX 



° cx dX ' 



(2.48) 



From Eq.(2.35) the solute molar fraction behaves as X = 

n 2 /n! = C(C + Citp)/n c0 - n 2c ^/n 2 cQ H at T = T c0 . 

For small T — T c0 and tj) it is expressed as 

= 1 + A m (T - T c0 ) - (K Kl ./n 2 c0 T c0 )yJ + • ■ ■ , (2.49) 

where A m is a constant. In two-phase coexistence this 
equation yields 



AX/Av = (K Kr /T c0 )X c 



(2.50) 



From Eq.(2.22) this is the near-critical expression of 
C(<9p/<9C)cx,t/T c0 = (dp/dA) c ^ T in the dilute limit. 

In Table 1, we show experimental data of T~^dT c j ' dX , 
(ricoTco^dpc/dX, (n^Tco)" 1 K Kr , and dp c /dX/K Kr for 
dilute mixtures near the solvent criticality, where the sol- 
vent is COz^ or H 2 02^. For C0 2 we have T c0 = 304K, 
n c0 T c0 = 26.1MPa, and {dp/ dT) cx /n c0 = 1.97, while for 
H 2 we have T c0 = 647.01K, n c0 T c0 = 96.0MPa, and 
(dp/dT) c Jn c0 = 1.81. Thus dTjdX, dpJdX, and K Kl . 
can be both positive and negative depending on the spe- 
cific details of the two components. These quantities are 



very small for H 2 0-D 2 mixtures-*', where the two com- 
ponent are very alike. If the solute is H 2 and the solvent 
is D 2 0, their signs are simply reversed with their absolute 
values nearly unchanged. 

In two-phase coexistence with general compositions, 
the present author introduced the parameter^ii, 



AX 

: An 



1 



dp 
OA 



(2.51) 



cx,T 



where n c = n lc + n 2c is the critical density and An = 
ri£ — n g . The critical line under consideration is that 
of the gas-liquid criticality for \e az \ ^ 1 and is that of 
the consolute criticality for \e az \ ^ 1- In the dilute limit 
X -> 0, we have e az S -(K K r/n c0 T c o)X . For 3 He- 4 He 
mixturesiii, the relation e az = —^X(l — X) roughly 
holds along the critical line, where X is the 3 He molar 
fraction. Thus / n c oT c $ is 1/3 with 3 He being a solute 
and is —1/3 with 4 He being a solute. Thus 3 He- 4 He 
mixtures are nearly azeotropic at any X (even away from 
the critical line). The resultant crossover effects have 
been observed in near-critical 3 He- 4 He mixtures in statics 
and dynamics^. 

On approaching the critical point, the thermal fluc- 
tuation of ip is enhanced with its variance propor- 
tional to the compressibility Kta — {dn/dp)T&/n as 
in Eq.(B5) in the appendix. As shown in Eq.(2.49) or in 
Eq.(B12), the thermal fluctuation of the molar fraction 
contains the growing part — (KKi/n 2 T c o)X4^^- From 
Eqs.(B5), (B8), and (B12) the concentration susceptibil- 
ity (8X/dA) p T behaves near the criticality as 

T {jk). - X + x2 ( K cr/ n coT c0 )K TA . (2.52) 

The first term is the low density limit (see Eq.(B8)). The 
second is the singular contribution stemming from the 
solute-solvent interaction. We may set (dX / 8A) p t = 
X and replace the mixture compressibility Kta by the 
pure-fluid compressibility when 



XK 2 CI K T /n c0 T c0 « 1. 



(2.53) 



This condition has been assumed in the definition of the 
Krichevskii parameter (see the appendix). 

3. Critical behavior of surface tension 

Using the Landau expansion of fo in Eq.(2.33) we next 
examine the mean-field critical behavior in two-phase co- 
existence, where the average order parameter values in 
the two phases are ip = ±ip e with 

vj e = L4o(T c0 -T)/B ] 1/2 . (2.54) 

The surface tension of the pure fluid 70 is written as 



7o = -{T c0 -T)A^lt 



(2.55) 
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TABLE I: T' c /T c0l p'Jn c0 T c0 , K Kl /n c0 T c0 , and p'JK Kl for 
CO2+ solute and for H2O + solute near the solvent critical 
point, where T' c = dT c /dX and p' c = dp c /dX. The last quan- 
tity is related to the temperature-derivative of the surface 
tension in Eq.(2.60). Data are taken from Refs.— >2S»2£. 



Solvent 


Solute 


T'JT c0 


p' c /ri c oT c o 


Kkt /TlcoTcO 


p'JK Kv 


C0 2 


Neon 


-0.0517 


0.919 


1.02 


0.900 


CO2 


Argon 


-0.192 


0.553 


0.936 


0.591 


CO2 


Ethanol 


0.539 


0.694 


-0.380 


-1.81 


C0 2 


Pentanol 


2.20 


1.96 


-2.42 


-0.809 


C0 2 


Ethane 


-0.182 


-0.187 


0.175 


-1.07 


H 2 


Toluene 


-1.32 


-0.948 


1.434 


-0.661 


H 2 


D 2 


-0.0050 


-0.0041 


0.0050 


-1.21 



The interface profile is expressed as ip(z) = ip e tanh(z/2£) 
along the surface normal, where £ is the correlation 
length in two-phase coexistence expressed as 

£ = (D n /2A )^ 2 (1 - T/Teo)- 1 / 2 . (2.56) 

Thus 70 oc (1 — T/T c o) 3 ' 2 , as originally derived by van 
der Waals^. 

It is easy to calculate the surface adsorption T in 
Eq.(2.28). Use of the expansion (2.35) gives T = 
<;C 2 fdz[i/j(z) 2 ~^ 2 e }. Thus, 



AT 

-2C 2 ^C = -2^4>V&, 



T, 



<::0 



(2.57) 



so T oc C(l-T/T c0 ) 1/2 . Because dj /dT = -3j /2(T c0 - 
T) = -2A ip 2 £ from Eq.(2.55), we find 



djo AT C 
dT T, 



<::0 



(2.58) 



If we write 70 = A s (l — T/T c0 ) 3 / 2 with A s being a con- 
stant, the surface tension of dilute mixtures 7 = 70 — TT 
is expressed as 



7 = A s T- 3/2 [T c (C)-T] 



3/2 



(2.59) 



to first order in £. That is, the solute effect on 7 is only 
to shift T c0 to T C (C) = T c0 + AT C . From Eqs.(2.31) and 
(2.58) we may express (dj/dT) cx in terms of AT C . Fur- 
ther using Eq.(2.42) it assumes a simpler form in terms 
of Ap c or dp c /dX as 



37 
dT 



1 djo Av Ap 

It ~ ~axt7 



r0 



1 dp c 
K Kl dX' 



(2.60) 



which tends to a well-defined limit at the solvent critical- 
ity. See the last column of Table 1 for the above ratio. 
It is negative if Kk t and dp c /dX have different signs. 



III. VAN DER WAALS THEORY OF MIXTURES 

A. Dilute mixtures 

The van der Waals theory of one-component fluids^ 
was extended to binary mixtures by van der Waals and 



Korteweg 3 ^. For binary mixtures the Hclmholtz free 
energy density / = f{n 1 ,n 2 ,T) is given by 



f 









-1 







ij 



(3.1) 



where Xi — (2Tr/miT) 1 ^ 2 h are the de Broglie lengthse 
with mi and m 2 being the molecular masses and h being 
the Planck constant. The cj) = fio«i+f20"2 is the volume 
fraction of the hard-core region with viq and v 2 q being 
the molecular volumes. The coefficients Wij represent the 
strength of the van der Waals attractive interaction be- 
tween ij pairs. However, more elaborate thermodynamic 
models have been used to predict the surface tension of 
real binary mixture a 18 i 19 ' 20 . 

In the pure fluid limit (n 2 = 0), the free energy density 
and the chemical potentials are given by 



fo(n,T) = Tnln 
Ho(n,T) = Tin 



1 - , 



— Tn — wiin 2 , (3.2) 



1 



2wnn, (3.3) 



where we set n = n\ and 4> = wqi^i- Hereafter 



e = v w wn 



(3.4) 



is the attractive energy among the molecules of the first 
component. In the pure fluid, the critical temperature, 
pressure, and density are written as 



Tc0 ~ 27' Pc0 ~ 27 Vl ° 



n c o = ^v 10 \ (3.5) 



See the upper plate of Fig.l for the liquid and gas den- 
sities in the van der Waals model. Far below the critical 
temperature in two-phase coexistence, the gas density 
n g o becomes very small compared to the liquid density 
riiQ. In fact, if p° cx = Tn g Q <C TngQ. the van der Waals 
theory yields 



<f> B /<f>t = (e^/T)e- £ ^( 2 -^)/ T , 



(3.6) 
(3.7) 



where <f>£ = vio««) is obtained from p° cx = and 4> g = 
v w n g0 from = Tln(n g0 A?). 

The quantity ip in Eq.(2.2) becomes 

<p = ^ - ln(l - 0) - ^w<t>, (3.8) 

in terms of (p = vqtii. Here two dimensionless parame- 
ters, the volume ratio and the potential ratio, are intro- 
duced as 



VQ2 
Vol ' 



Wl2 



(3.9) 



which characterize the physical properties of the second 
component. If n(z) — n\(z) is the density profile of the 
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pure fluid across an interface, the density ni is expressed 
as in Eq.(2.9). With the aid of Eqs. (2.7) and (3.3) we 
rewrite (p as 



1 o 

9 = yMc 



ln(nA?) + Dun" 



+ - 



—— ( w 



(3.10) 



in terms of cf> = vion(z). From Eq.(2.9) the space- 
dependent molar fraction X(z) — ni{z)/ni(z) becomes 



X = ( exp 



i_l0 + ^(«,-l)0 + D'n" 



(3.H) 



where C = Afe^-/ T C- = (m 2 /mi) 3 / 2 e(' ia -''°W T and 
D' = D12 — Dn. Notice that X =const. or JC = 1 for 
r = 1, w = 1, and D12 = -D11, where the two components 
have the same physical properties. 



B. Two-phase coexistence 

From Eq.(3.11) the logarithm of the partition coeffi- 
cient K, in Eq.(2.24) is expressed as 



ln/C 



(r-l)A^) 2e 
(l-&)(l-<£ fl ) ~ T 



(w - l)A(f>, (3.12) 



where <j> e = v 01 n m , 4> g = v in g0 , and A(/> = (f>e - (f> g = 
vioAn. In the lower plate of Fig.l, /C vs T/T c q is shown 
for typical four cases. Remarkably, the azeotropy (/C = 1) 
is attained in the dilute limit on the following line in the 
r-w plane, 




0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 

T/T c0 

10r r- 




0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 

T/T c0 

FIG. 1: Upper plate: A Z (T) in Eq.(3.13), n l0 /n c0 , and 
n g o/n C Q vs T/T c o in the van der Waals theory. Lower plate: 
K, — X g /X£ vs T/Tco in dilute mixtures on a semi-logarithmic 
scale, where (r,w) = (0.8,0.8) for (a), (0.8,1.0) for (b), 
(1.5, 1.0) for (c), and (1.5, 1.4) for (d). Two parameters r 
and w are defined in Eq.(3.9). 



r-l = A z (T)(w-l), (3.13) 

where the coefficient A z (T) is determined by the solvent 
properties only as 

A Z {T)= 2(1 -ct> t ){l-^ g )ejT. (3.14) 

See the upper plate of Fig.l for A Z (T) vs T/T c0 . Here 
A z — > 3 as T — > T c q, while for <f> g 1 we find A z = 2/<f>i 
from Eqs. (3.6) and (3.7) and 

ln/C [r- 1-A z (w -l)]$e/T. (3.15) 

Using kn in Eq.(2.26), Sengers el alX examined kn/fo 
where /o = exp(/io/T) is the fugacity of the pure fluid. 
In the van der Waals theory it is of the form, 

k H 3 mi r- 1 2e 

In— = -ln — + - — <t> t --= (w-l)fa. (3.16) 

Jo 2 m 2 1- 4>£ T 

In Fig. 2, we display profiles of the molar fraction X(z) 
divided by the molar fraction in the gas region X g at 
T = 0.9T C , where A z = 2.88. In the left panel, we set 
r = 0.8 and vary w as 1.3, 1.1,0.8, and 0.5. In the right 
panel, at r = 1.5, we have w = 1.5, 1.2, 0.8, and 0.4. Thus 



Xi/X g = K 1 increases with decreasing r and/or with 
increasing w. In our numerical analysis, we set D u = 
D12 = 10a 5 with a = v\(^ . See the next subsection for 
justification of this choice of Du. It is worth noting that 
Sahimi and Taylor calculated the density profiles of two 
components around an interface^. 

C. Surface tension 

In the upper plate of Fig. 3, we show 70 and 7o/t 3 / 2 vs 
T/T c0 for the pure fluid, where we used the formula (A5) 
in Apenndix A with Dn = 10a 5 . The relation 70 oc r 3 / 2 
nicely holds over a wide range of T/T c q. Remarkably, ex- 
perimental data of the surface tension of water can also 
be nicely fitted to the formula 7 oxp = A exp (l — T/T c0 ) 3 / 2 
over a wide temperature range except close to the criti- 
cality (in the range 1— T/T c0 > 0.1)2i. As in our previous 
worb 2 ^, we have determined Dn such that our numerical 
7 and the experimental 7 cxp for water reasonably agree 
except close to the criticality. In fact, at T/T c = 0.675, 
our 7 is 42.5 dyn/cm if we set Dn = 10a 5 , a = 3A, and 
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r =0.8 r = 1 .5 




z/a z/a 



FIG. 2: Profiles of the normalized molar fraction X(z)/X g 
of the solute at T = 0.9T C for r = 0.8 and w = 1.3, 1.1,0.8, 
and 0.5 (left) and for r = 1.5 and w = 1.5,1.2,0.8, and 0.4 
(right). It tends to K,^ 1 on the left hand side. Shown also is 
the normalized density von(z) of the pure fluid composed of 
the first component (broken line). Space is measured in units 

1/3 

of the molecular size a = v-, A . 




°' 2 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 



T/T c0 

FIG. 3: Upper plate: Surface tension 70 multiplied by 10 and 
70/V 3 / 2 vs T/T c o for the pure fluid. Lower plate: A'y/Xg = 
-TYjXg vs T/T c o in dilute mixtures, where (r, w) = (0.8, 0.8) 
for (a), (0.8,1.0) for (b), (1.5,1.0) for (c), and (1.5,1.4) for 
(d). The 70 and A7 are measured in units of e/a 2 . 



r = 0.8 r = 1.5 




260 270 280 290 300 ' 2 60 270 280 290 300 

z/a z/a 



FIG. 4: Excess solute density Ari2(z) in Eq.(3.17) divided 
by the molar fraction in the gas X g in units of v± . 




0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 



T/T c0 

FIG. 5: Ratio of the temperature derivatives 7^ x = 
((97/(9T) cx ,p and 70 = dr/o/dT vs T/T c o in dilute mixtures, 
where (r,w) = (0.8,0.8) for (a), (0.8,1.0) for (b), (1.5,1.0) 
for (c), and (1.5, 1.4) for (d). It changes its sign for (d). 



T c q = 647. IK, while the experimental value of water is 
44.6 dyn/cm. 

In the lower plate of Fig. 3, we display the surface ten- 
sion change A7 = — TV divided by X g vs T/T c q for four 
sets of (r, w). From Eq.(2.27) T is the space integral of 
the excess solute density Ari2(z) expressed in terms of 
the density n(z) of the reference pure fluid, 



An 2 (z) = n 2 (z) - n 2 i (n(z) - n ea ). (3.17) 

In Fig. 4, we plot n 2 (z) for r = 0.8 (left) and r = 1.5 
(right) for various w. With increasing u>, T becomes neg- 
ative and its magnitude increases strongly. In Fig. 5, we 
display the ratio (dp c /dX)/K Kl = {d-//dT) c ^ p /{d-//dT) 
calculated from Eq.(2.31) as a function of T/T c q for four 
sets of (r, w) . It even changes its sign from positive to 
negative with increasing T for (r, w) = (1.5, 1.4). 



J_ d_P c 
n c0 T c0 d X 



0.5 1 1.5 2 2.5 3 

W 




FIG. 6: T^dTc/dX ( upper plate) and (n c oT c o) 1 dp c /dX 
(lower plate) as functions of r and w in Eq.(3.9) near the 
solvent critical point. 



D. Near-critical behavior 

The Landau expansion of /o with respect to ip = n—n C Q 
is given in Eq.(2.32). For the van der Waals model the 
coefficients are given by 



27 243 o 

M = — v w , B Q = — T c0 v 10 . 



(3.18) 



In the pure fluid, the liquid and gas densities are n^o 
n c o + ijj e and n 9 o = n C Q — ip e , where Eq.(2.49) gives 



%jj e = 2n cQ T 



1/2 



(3.19) 



Here r = 1 — T/T c o is the reduced temperature (positive 
below the critical temperature). See the upper plate of 
Fig.l for neo and n g o. Then, 

An = 2ip e , Av = 18vf ip e , As = 9v w ^ e . (3.20) 

These differences are of order t 1 / 2 . In particular, As = 
6T 1 / 2 . The latter two relations are consistent with 
the Clausius-Clapeyron relation As/Av — {dp/dT) cyi = 
l/2uio along the coexistence curve. In addition, the cor- 
relation length £ in Eq.(2.51) becomes £ = 0.86ar -1 / 2 in 
our numerical analysis with Dn = 10a 5 . 

Using if in Eq.(3.6) we perform the Taylor expansion 
e-f = C + Clip + C 2 ^ 2 /2 + • ■ • as in Eq.(2.36). In terms 



FIG. 7: Curves of AT C = 0, Ap c = 0, and AX = on the 
plane of (r, w) near the solvent critical point. 



of r and w in Eq.(3.9) the coefficients are expressed as 
Co - ^exp(-r/2 + 9w/4), 

Ci = ^v w (-3r - 2 + 9w)C , 
81 

C 2 = — v 2 10 {{r - 3w) 2 - 4w}C . (3.21) 
The critical solute density and concentration are 

n 2c = CoC, X c = 3vioCoC. (3.22) 

The solute density difference in Eq.(2.18) and the com- 
position difference in Eq.(2.19) are expressed as 



A?i 2 = (3r - 9w + 2)X c n c0 T 1/2 , 
AX = 3(r-3w + 2)X c r 1/2 . 



(3.23) 



The Krichevskii parameter in Eq.(2.47) is given by 



K Kr /n cQ T cQ = -(r - 3w + 2), 



(3.24) 



which was already derived by Petsche and Debenedetti^. 
See Table 1 for experimental values of the above quantity. 
In accord with these results, K behaves as 



K= l + 3(r-3w + 2)T 1/2 + ---. 



(3.25) 



while ln(fcff// ) = const. + 3(r-3w + 2)r 1 / 2 /2H from 

Eq.(3.16). Sengers et alA found that data of T\n(kH / fo) 
can well be fitted to the form C + B(ni — n c o) near the 
critical point for a number of solutes in H 2 0. 

From Eqs. (2.43) and (2.44) the derivatives dT c /dX and 
dp c /dX along the critical line are written as 



1 dT c _ 1 
7k dX ~ 4 

1 dp c 
Pco dX 



(r - 3w) 2 - to, (3.26) 
(r-3tu)(r-3to-|-2)-4(u>-l). (3.27) 
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In Fig. 6, we show dT c /dX and dp c /dX in the r-w plane. 
In Fig. 7, we show the curves of dT c /dX = 0, dp c /dX = 0, 
and the azeotropic line AX = 0. Thus, AX, dT c /dX, 
and dp c / dX can be both positive and negative depending 
on r and w. 

From Eqs.(2.58) and (3.26) the surface adsorption T is 
written as 

T=\[(r--iwf-MX c d ^. (3.28) 

From Eq.(2.60) we calculate the temperature-derivative 
of 7 on the coexistence surface, 



d-f 
df 



cx,p 



3w 



2(w - 1) 



r - 3w + 2 



Thus the above derivative can be both negative and pos- 
itive and can even diverge to ±oo on the azeotropic line 
r - 3w + 2 = 0. 



IV. SUMMARY 

In summary, a Ginzburg-Landau theory has been 
presented for dilute binary mixtures, where the solute- 
solvent interaction is relevant but the solute-solute 
interaction is negligible. A parameter £ proportional 
to the solute fugacity has been introduced in Eq.(2.8). 
Up to first order in £, all the physical quantities of 
binary mixtures can easily be calculated in terms of 
the properties of the one-component fluid and the 
solute-solvent interaction parameters. In more detail, 
our main results are as follows 

(i) The coexistence surface has been given by Eqs.(2.16) 
and (2.23) or by Eq.(2.22). Henry's constants have been 
introduced in Eqs.(2.24)-(2.26). 

(ii) The Gibbs formula for the surface tension change 
A7 in Eq.(2.27) has been derived in Appendix A. The 
surface tension derivative (dj/dT) pc ^ with respect to T 
at fixed p has been obtained in Eq.(2.31). Interestingly, 
it consists of two terms both being independent of (. 

(iii) The critical temperature shift AT C is given 
in Eq.(2.41) and the critical pressure shift Ap c in 
Eq.(2.42). 

(iv) The Krichevskii parameter Kk t has been given 
in Eqs.(2.47) and (2.48). The normalized parameter 
Kj<irX c /n c oT c o represents the size of the critical concen- 
tration fluctuations as in Eq.(2.49), leading to Eq.(2.52). 

(v) The surface adsorption T has been realted to AT C 
as in Eq.(2,58) and {d"f/dT) c ^ p to Ap c as in Eq.(2,66) 
near the criticality. The solute effect on the near-critical 
surface tension is simply to shift the critical temperature 
T c0 by AT C as in Eq.(2.59). 

(vi) Experimental data of dT c /dX, dp c /dX, and -Kxr 
have been given in scaled forms in Table 1, which shows 
that they can be both positive and negative. 

(vii) The van der Waals model of binary mixtures 
has given simple expressions for all the theoretical 



expressions in Section II, as illustrated in the figures. 
The solute-solvent interaction is described in terms of 
the size ratio r and the potential ratio w in Eq.(3.9). 

(viii) The profiles of the solute density and its excess 
near an interface have been numerically calculated as in 
Fig. 2 and 4. The negative adsorption becomes marked 
for large w. 

(ix) The near-critical behavior in the van der Waals 
model is very simple in the mean-field theory. In 
terms of r and w we have calculated JC, Kk t , dT c /dX, 
dp c /dX, T, and (dj/dT) c ^ p . In the r-w plane , we have 
plotted dT c /dX and dp c /dX in Fig. 6 and the curves of 
dT c /dX = 0, dp c /dX = 0, and K = 1 in Fig. 7. 

Finally, we propose measurements of the surface ten- 
sion as a function of the temperature in the isobaric con- 
dition for various solutes in water or in C02- From our 
theory, the derivative (dj / '9T) cx . p becomes independent 
of the solute density in the dilute limit and can be both 
negative and positive, being delicately dependent on the 
size ratio r and the potential ratio w in the van der Waals 
theory as in Eq.(3.29). It is a relevant parameter deter- 
mining the Marangoni flow around a bubble moving in 
heat flow in binary mixtures, as will be reported shortly. 

Appendix A: Calculation of surface tension 

Here the surface tension 7 of binary mixtures is ex- 
amined from Eq.(2.1). The grand potential density of 
mixtures is given by 



T 



/ / 

Hill 3 , 



(Al) 



1,3 



where n- = dru/dz and fii take the values in two-phase 
coexistence. All the quantities change along the z axis. 
The space integral of uo gives the grand potential £1 in 
Eq.(2.10). Then ui tends to — p cx far from the interface 
z —> ±00 and the surface tension is expressed as 



7 



/ 



dz[uj(z) +p c 



(A2) 



Differentiation of ui(z) in Eq.(Al) with respect to z 
yields dw/dz = 2/ V (; /■>,,/'!"'' from Eq.(2.6), where 



d rij/dz . Therefore, 



n'i n 'j 



2(/ - ^fiiUi) +p c 



(A3) 



Then 7 in Eq.(A2) may also be expressed as 7 = 

/ dzTJ^ij D ij n i n j = 2 f dz (f - Si £» n * +Pcx)- 

Next 7 is expanded with respect to £ in the dilute 
case. As in the derivation of Eq.(2.11), elimination of fa 
in Eq.(Al) gives 



7 = / dz 



/o(m)- 



■fiWi +p C y i + —D 11 ni 



12 



■Tno 



(A4) 
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where fix — + 6^i lcX7 n 2 is given by the second line 
of Eq.(2.9), and the integrand vanishes as z — > ±00. Let 
n = n(z) be the density of the reference pure fluid or 
n(z) = lim^o rii(z). Then n{z) — > ri£Q (n g o) as z — > —00 
(00) and we have the interface equation (2.7). As £ — > 
the surface tension of the pure fluid is obtained as 



To 



/o(n) - M° x n- 



(A5) 



From Eqs.(A4) and (A5) the surface tension change 
A7 = 7 — 70 for small Q is expanded with respect to 
the deviation 5ni(z) — ri\(z) — n(z) as 



A7 



dz 



(fo(n) 



/4 



TD u n")Sm 



-Snivji + Sp cx - Tn 2 



+ 



(A6) 



to first order in £. Here the first term in the brackets 
vanishes from Eq.(2.7). Further use of Eqs.(2.16) and 
(2.17) yields the Gibbs relation in Eq.(2.27). 

Appendix B: Correlation- function expressions 

We examine the correlation-function expressions for 
thermodynamic derivatives such as K^r in Eq,(2.47) and 
(dX/dA) p T in Eq.(2.52) in the framework in the book of 
the present author—. Equivalent relations for Kkt were 
already used in the literature ^ 14 ' 15 . 

The microscopic particle densities are written as 



eej 



rt), 



(Bl) 



where the summation is over the particles of the species 



j(= 1,2) at position re. Then rij 
denotes the equilibrium average, 
functions are written as 



= (hj), where (• • •) 
The pair correlation 



(6ni(r)6nj(Q)) = niSijS(r) + riinjgij(r), (B2) 

where 8hi(r) = hi(r) — rij (i = 1, 2) are the density devi- 
ations and gij(r) = 1,2) are the radial distribution 
functions tending to zero for large separtion r. It is con- 
venient to introduce the concentration variable X(r) and 
the number density variable h(r) by 



X = X 



jn 2 



n 2 



(B3) 



where (h) = n = n% + n 2 and (X) = X = n 2 /n. We de- 
fine a fluctuation variance for any space-dependent vari- 
ables A(r) and B{r) by 

(A:B) = J dr((A(r) - (A))(B(r) (B))). (B4) 

The variances among n and X may be expressed in terms 
of the thermodynamic derivatives, 

(dn\ ~ - T(dX\ 



TA 



(B5) 



where n and X are treated as functions of the field vari- 
ables T, p, and A = /i 2 — f^i in the derivatives. These 
variances are linear combinations of the variances among 
the densities, which are written as 



(B6) 



from Eq.(B2). On the other hand, the compressibility at 
constant X is written as 



Ktx 



1 fdn 



n\dp 



TX 



n 2 T 



{h : n) — 



(h : Xf 



(X:X) 



(B7) 



Near the mixture criticality, the ratio K^x /Kta behaves 
as X/(X : X) = nX/[(dX/dA) Tp T] (see Eq.(2.52)). All 
the variances in Eqs.(B5) and (B6) diverge strongly at 
the mixture criticality except for special cases such as 
the critical azeotropy. In the low density limit X — > 
under Eq.(2.44), Eqs.(B3) and (B6) give 



(X : X) Si (h 2 : n 2 )/n 2 S X/n. 



(B8) 



We also need to assume (X : n) oc X for the existence of 
the Krichevskii parameter (see Eqs.(B9) and (BIO)). 

We next examine the thermodynamic derivative 
(dp/dX) n T = —(dn/dX)pT/nKTx- Its correlation- 
function expression reads 



dp 
dX 



-nT(h : X) 



(n : n)(X : X) - (n : A) 2 



(B9) 



In the low density limit we use Eq.(B8) and replace the 
denominator of Eq.(B9) by (h : n)X/n to find 



hm-(H 

x^o nT\dX /nT 



lim 



n{h : X) 



x^o X{n : n) 
= 1-mCg, (BIO) 



where the second line follows from Eq.(B3). We define 



Co 
1 



lim (n 2 : ni)/n 2 (n 1 : hi), 

U2- >0 



(Bll) 



which coincides with the space-integral of the direct cor- 
relation function C\ 2 (r) in the dilute limft£ ' 14 i 15 . Here 
we define Gij{r) in dimensionless formsi Thus the 
Krichevskii parameter i^Kr in Eq.(2.47) is the value of 
n\T{l — niC^) at the solvent criticality. This expres- 
sion has been used to estimate for given molecular 
interaction parameter o 6 ' 14 ' 15 . From Eqs. (BIO) and (Bll) 
the singular parts of X and n 2 are 



(A) sing = (C~ - l/n^XSnt 
(n 2 ) S m g = C^n 2 Shi, 



(B12) 



near the mixture criticality. Here we have calculated pro- 
jected parts of X and h 2 onto the critical fluctuation 
8fii = hi — rii- Equation (2.52) is then obtained with 
the aid of Eq.(B8). 
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